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We discuss nuclear modification of fragmentation functions in the context of the so- 
called "rescaling models". These models implement partial deconfinement inside nuclei 
by modifying the fragmentation functions perturbatively. We apply these models to the 
analysis of nuclear hadron production in deep inelastic scattering processes at the HER- 
MES and EMC experiments. 

In deep inelastic scattering a projectile lepton I emits a virtual photon 7*, which scatters 
on a quark q from the target T, in our case a deuteron Dora heavier nucleus with atomic 
number A. The struck quark fragments into the observed hadron h, see Fig. [TJ We use 
kinematic variables as summarized in the same figure: x is the Bjorken's scaling variable, 
v is the energy of the virtual photon in the target rest frame, and z is the fraction of the 
energy transferred to the produced hadron. 

The experimental data on nuclear effects in hadron production are usually presented in 
terms of the multiplicity ratio 

where N l A is the number of outgoing leptons in DIS processes on a nuclear target of atomic 
number A and dN\/dz is the ^-distribution of produced hadrons in the same processes; 
the subscript D refers to the same quantities when the target is a deuteron. A similar 
definition for the multiplicity ratio as a function of v may be written. In absence of nuclear 
effects the ratio R h M would be equal to 1 . The experimental observation that R h M 7^ 1 [l|-|3| 
has been explained theoretically in many ways: as an effect of nuclear absorption of the 
produced hadrons ||||, in a gluon-bremsstrahlung model for leading hadron production 
as an higher- twist effect J7|. In this short note we will not discuss these models in 
detail (see Ref. ||) and, instead, will concentrate on the so-called rescaling models pHII 
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Figure 1. Kinematic variables for DIS scattering; four-moments are 
indicated in parentheses. 



The starting point of rescaling models is to assume a change in the confinement scale 
Xa in nuclei, compared to the confinement scale Ao in free nucleons: 

A a > A • 

This assumed partial deconfmement in nuclei affects both the parton distribution functions 
(PDF) and the fragmentation functions (FF). For what concerns PDF's consider 



a valence quark which carries a momentum Qo when it is confined on a scale Ao- If the 
scale changes to Xa it carries a corresponding momentum Qa- Since there is no other 
dimensionful scale, the product QX must remain constant, so that 

QoAo = QaXa ■ (2) 

Therefore, if we take Qo to be the initial scale for the DGLAP evolution of distribution 
functions, we may set 

qf (x,Qa = T^Qo) = Qf(x,Qo) , (3) 
where qf is the distribution function of a quark of flavour / in a free nucleon and is the 



same quantity when the nucleon is inside a nucleus. For FF's |IT] a similar argument holds 
because constituent quark and effective hadron masses are sensitive to the confinement 
scale. Therefore if we take Qo to be the physical threshold for hadron production we may 

set 

D t f A (x, Qa = -p" Qo) = Df(x, Q ) , (4) 

AA 

where Df is the fragmentation function of a quark of flavour / into a hadron h and D h f A is 
the nuclear modified fragmentation function. To extend Eqs. (|3|) and (^) to an arbitrary 
scale Q we apply pQCD evolution. The nuclear structure and fragmentation functions 
evolve over larger interval in momentum compared to the corresponding functions at the 
same scale Q, since the starting scale is smaller, see Eq. (0). The final result [|g|— fTTH is 

qf(x,Q)=q f (x,U(Q)Q) (5) 

D h f lA (z,Q)=D h f (x,UQ)Q) , (6) 

where the scale factor £,a{Q) is greater than one. These two equations are the main tool 
which we will exploit. There are two models for the computation of the scale factor: 
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Figure 2. Comparison of partial and maximal dconfinement models at EMC with a copper target. 
a) ^-distribution with rescaling only and b) with rescaling and absorption; c) z-distribution with 
rescaling and absorption. Data taken from Ref. 0] . 

a) the maximal deconfinement model (NP) || assumes the onset of "colour conductivity" 
in nuclei and takes A^ = Ra, where Ra is the nuclear radius, so that 

£,a(Q) = Ra/R p ; 

b) the partial deconfinement model (CJRR) Jl0[ assumes the deconfinement scale A^ to 



be proportional to the degree of overlap of the nucleons inside the given nucleus, and 

. 1 «s(Qq) 

UQ) = (Aa/a ) 3 - (Q) • 

Note that the maximal deconfinement model assumes a much larger scale factor, which 
results in a larger nuclear modification of hadron production at high z. 

We calculate the multiplicity ratio ([I]) using the rescaled PDF ([]) and the rescaled FF 
@ in the leading order pQCD computation of N A and dN A /dz: 



1 dN\ 
~W A ~dz~ 



— / dxduY,e 2 f q f (x,UQ)£-r-D h f (z,UQ) 

u Jexp. cuts f tlxL LLiy 



(7) 



a 



7*A 



dx dv e 2 f qf(x, £a Q) 

exp. cuts j 



do-"<* q 
dxdv 



where e/ is the electric charge of a quark of flavour /, da 1 q /dx dv is the differential cross- 
section for a 7*g scattering computed in pQCD at leading order, and Q 2 = 2Mxv, with 
M the nucleon mass. In the numerical computations we used GRV98 parton distribution 
funtions \12\ and Kretzer's parametrization of FF's at leading order [[0| . 

The results for the ^-dependence of the multiplicity ratio for charged hadrons are shown 
in Fig. ||a and compared to EMC data ||. The maximal deconfinement model underes- 
timates the data over nearly the whole range. On the contrary the partial deconfinement 
model fit the data at high is, but overestimates them at smaller is. Because of Lorentz 
dilatation, at high is the hadron is formed mainly outside the nucleus and it is affected by 
rescaling effects only On the contrary, at small is the hadron is formed inside the nucleus 
and starts interacting with its nucleons, with some probability of being absorbed. 

To take into account nuclear absorption, we follow the analysis of Ref. |4]|J: the quark 
propagates in the nucleus with a cross-section a q for inelastic scattering on the nucleons. 
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Figure 3. Comparison of absorption effects in z-distributions with rescaling only (dashed lines) and 
with rescaling and absorption (solid lines). From left to right: EMC with copper target, HERMES 
with krypton and HERMES with nitrogen targets. Data taken from Refs. Note that the scale 

in Fig. ||a is different from Fig. ^c. 



Subsequently, it creates a pre-hadronic states, which has a cross-section cr*, and finally 
the observed hadron is formed, which has a cross-section <jh- Following the experimental 
indications of EMC and HERMES, we take o~ q = and <J* = = 20 mbarn for charged 
hadrons. With the latter assumption we need to consider effectively only the formation 
time t f for the pre-hadronic state 0: 



tf 



— hi(z 2 ) — 1 + z 2 \ zv 



1 — Z 2 J K 



where k is a parameter which we fix to k = 0.25 GeV/fm. Note that rp— >--(l — z)v as 
z—*l, giving the formation time suggested by the gluon bremsstrahlung model of Ref. |J. 
Finally, nuclear absorption effects are included in the computations by multiplying the 
integrand in (|^) by the nuclear absorption factor A/a, which represent the fraction of 
hadrons which escape from the nucleus: 

POO 

2i. / » _ / l ..\ ro ( u ..mA-i 



■A/a = / d 2 b / dyp A (b,y)[S A (b,y)} 

J J — oo 

where p A {b,y) is the nuclear density normalized to 1 at transverse and longitudinal coor- 
dinates (b,y), and SA(b,y) is the survival probability of a hadron produced in (b,y): 



S A (b,y) = l-a h ^ dy'p A (b,y')(l-e-^-yy^ . 
For the deuteron we used as a density the sum of the Reid's soft-core S- and D-wave 



functions squared |14|| . For heavier nuclei we used a Woods-Saxon density with radius 
R A = 1.12A 1 / 3 - 0.86A- 1 / 3 fm. 

The v- and z-dependence of the multiplicity ratio at EMC and HERMES after the 
inclusion of nuclear absorption are shown in Fig. [2]b and Fig. |2|c. It is clearly seen that 
while the CJRR model gives a nice description of the data, the NP model is ruled out. 

The average values of (v) available at HERMES are smaller than at EMC and the 
absorption effects larger. In Fig. |3] we compare the z-distributions in the partial decon- 
finement models at EMC and HERMES before and after inclusion of absorption. As 
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expected, at EMC with a 63 Cu target nuclear absorption is marginal and rescaling alone 
gives a satisfactory description of the data. At HERMES with a 84 Kr target, which is 
comparable to copper, both effects are larger and absorption is dominant, tending to 
mask rescaling effects. With a 14 N target, which is smaller than krypton, both effects are 
smaller. 

In summary, rescaling models (supplemented by nuclear absorption) are shown to be 
able to describe both HERMES and EMC data on the nuclear modification of hadron 
production in DIS processes. While the maximal deconfinement model is ruled out by the 
data - it assumes a too large deconfinement - the partial deconfinement model is shown 
to be a good one. Further precise data at moderate and high z/s and for light and heavy 
targets are needed to disentangle rescaling and formation time effects. 
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